Classification of time-reversal-invariant crystals with gauge structures

A peculiar feature of quantum states is that they may embody so-called projective representations of symmetries rather than ordinary representations. Projective representations of space groups-the defining symmetry of crystals-remain largely unexplored. Despite recent advances in artificial crystals, whose intrinsic gauge structures necessarily require a projective description, a unified theory is yet to be established. Here, we establish such a unified theory by exhaustively classifying and representing all 458 projective symmetry algebras of time-reversal-invariant crystals from 17 wallpaper groups in two dimensions-189 of which are algebraically non-equivalent. We discover three physical signatures resulting from projective symmetry algebras, including the shift of high-symmetry momenta, an enforced nontrivial Zak phase, and a spinless eight-fold nodal point. Our work offers a theoretical foundation for the field of artificial crystals and opens the door to a wealth of topological states and phenomena beyond the existing paradigms.

Symmetry groups and their representations are at the heart of physics. When going from classical to quantum physics, a classical symmetry group G becomes represented in the Hilbert space, where it makes no physical difference if all states are multiplied by a global phase. It follows that the representation allows an extra phase factor, i.e., for g 1 , g 2 ∈ G, their representations ρ(g 1 ) and ρ(g 2 ) may satisfy ρ(g 1 ) ρ(g 2 ) = ν(g 1 , g 2 )ρ(g 1 , g 2 ) with ν(g 1 , g 2 ) ∈ U(1). This is known as the projective representation of G, and the phase factors ν called the factor system for this representation. As a well-known example, classifying the projective representations of Poincaré group for elementary particles leads to the two types of particles, bosons and fermions, corresponding to two distinct factor systems 1 .
The defining symmetries for crystals are the space groups. What is the physical meaning of a projective representation in this context? Consider a spinless quantum particle on a two-dimensional (2D) lattice as shown in Fig. 1a. A projective representation for the lattice translations allows (hereafter, we use bold letters to denote the represented symmetry operators) L a L b L À1 a L À1 b = e iθ , from which one observes that the extra phase factor corresponds to a gauge flux through the lattice. This shows that projective representations of space groups are associated with gauge fluxes, and somewhat explains the previous ignorance of them in textbooks on solid-state physics 2 . Because most electronic crystals are free of gauge flux, one can show that their descriptions are restricted to ordinary representations. Nevertheless, it was recognized that rich gauge-flux configurations can emerge in certain strongly correlated spin systems, where projective representations of space groups are needed for their description [3][4][5][6][7][8] .
The rise of artificial crystals in recent years completely changes the situation. Artificial crystals cover a wide range of systems, such as acoustic, photonic, mechanical, circuit, and cold-atom systems [9][10][11][12][13][14][15][16][17][18][19] . Most artificial crystals intrinsically preserve time-reversal (T) symmetry, which allows fluxes 0 and π over the lattices. A salient feature is that these lattice gauge fluxes can be readily engineered. Recent works showed that these fluxes modify the physics in a fundamental way and projective representations are indispensable for understanding artificial crystals [20][21][22][23][24][25] . This urgently calls for a unified theory of projective representations of symmetries for T-invariant crystals, which constitutes the foundation of the whole field.
In this work, we develop such a theory and predict its distinguishing consequences. First, we characterize all possible projective symmetry algebras (PSAs) with time-reversal symmetry for any space group. This is demonstrated by 458 PSAs-189 of which are algebraically independent-for all 17 wallpaper groups in two dimensions. Then, we show all the 2D PSAs can be systematically realized by lattice models with appropriate gauge fluxes. Finally, we present three signature results of PSAs, including the shift of high-symmetry momenta, an enforced nontrivial Zak phase, and a spinless eightfold nodal point.

Projective symmetry algebras with time-reversal invariance
We start by presenting a general result that reduces the problem for systems with time reversal symmetry T. Let G be the space group, then the total symmetry group is G × Z T 2 , where Z T 2 = fE,Tg is the twoelement group generated by T. Mathematically, the classification of all possible factor systems for this group corresponds to the second group cohomology H 2,c ðG × Z T 2 ,Uð1ÞÞ 26 . We have proven that due to the anti-unitary character of T, H 2,c ðG × Z T 2 ,Uð1ÞÞ ffi H 2 ðG, Z 2 Þ × H 2 ðZ T 2 ,Z 2 Þ with Z 2 = f ± 1g (see Methods). Hence, the computation is simplified to deriving H 2 ðG, Z 2 Þ, since it is known that H 2 ðZ T 2 , Z 2 Þ ffi Z 2 = f ± 1g distinguishes integer and half integer spins, respectively. This means we only need to consider G with factors restricted to Z 2 = f ± 1g. Our discussion below focuses on spinless systems, which are pertinent to most artificial crystals. It can be directly extended to spinful systems, as we shall comment at the end.
The second group cohomology H 2 ðG, Z 2 Þ can be derived from the abstract group cohomology theory, e.g., from the twisted tensor product of the cochain complex of the translation subgroup and that of the point group 27 . Considering the 17 wallpaper groups for 2D, the results are listed in the second column of Table 1. Besides the classification, we also need to know the content of each H 2 ðG, Z 2 Þ, namely, the concrete algebraic relations satisfied by the symmetry operators, which are called the PSAs, because they are directly related to the physics of a system. We have worked out all PSAs in terms of generators of each group, as listed in the fourth column of Table 1. The technical details are given in Supplementary Note 2. The classification is complete, meaning that any T-invariant crystal system in two dimensions must belong to one of the PSAs listed here.
Interestingly, each H 2 ðG,Z 2 Þ is a product of Z 2 , i.e., H 2 ðG,Z 2 Þ ffi Z n 2 . Meanwhile, we find that the corresponding PSA can be captured by a complete set of nZ 2 -valued cohomology invariants, which are denoted by σ, α, β, η, and τ in Table 1. The specific meaning of these symbols will be explained in a while. Here, one can easily check that they are indeed cohomology invariants, by noting that they are unchanged when multiplying symmetry operators by arbitrary Z 2 phases ± 1. The ordinary representation just corresponds to the case with all invariants being + 1.
It should be noted that for each G the PSAs classified by H 2 ðG,Z 2 Þ have redundancies for their abstract algebraic structures, because often G has nontrivial automorphisms such that two nonequivalent factor systems lead to equivalent algebraic structures. We screen out all nonequivalent algebras, of which the numbers (N G ) are listed in the last column of Table 1. We find that there are 189 non-equivalent algebraic structures out of the 458 PSAs.
To illustrate our theory, we take group P2 as an example. The set of generators of P2 consists of two unit translations L a , L b and the rotation R by π (along an out-of-plane axis). Its group algebras are expressed in terms of the four combinations, R, L a R, L b R and L a L b R, each of which is squared to 1. According to Table 1, H 2 ðP2, Z 2 Þ ffi Z 4 2 , so there are four cohomology invariants α i (i = 1, 2, 3, 4), corresponding to the four PSA relations: Since any permutation of the four twofold rotations above gives an isomorphic PSA, there are only five equivalence classes of PSAs, and each class is specified by how many α^'s equal − 1.
Our result shows that the 17 wallpaper groups together with T symmetry can generate 458 × 2 PSAs (the additional factor 2 is from spin), which is much richer than the case of Poincaré group with only twofold classification. Physically, this is due to the reduced symmetry which allows more gauge flux configurations.

Flux realizations of projective symmetry algebras
After completing the classification, our next task is to develop a construction method to realize each of the PSAs. This is important for two purposes. First, it serves as a validity check for our results in Table 1 and a Gauge flux in a lattice requires the space group symmetries to be projectively represented. b Work flow diagram of this work. Given a space group G, the classification of projective representations are given by H 2 ðG,Z 2 Þ. These representations are captured by projective symmetry algebras (PSAs) with a complete set of cohomology invariants, from which we construct a canonical model that can realize all possible PSAs, and derive nontrivial physical consequences. c-h Illustrate the construction method for realizing the five basic classes of PSAs (cohomology invariants). Specifically, (c, d, e, h) are for nontrivial σ, α, β, τ, respectively, and (f and g) are for nontrivial η. In (g), the bond connecting the two black sites can be either present or absent, and the total flux through the vertical and horizontal rectangular plaquettes is required to be π. Here, the π flux in a plaquette is realized by a negative hopping amplitude (marked by red color) on an edge of that plaquette.
demonstrates that each PSA can indeed be realized in a physical system. Second, it provides guidance for the experimental realization of nontrivial PSAs in artificial crystals.
Our construction is via mapping each PSA in Table 1 to a specific gauge flux pattern. In this process, we distinguish five classes of cohomology invariants in PSAs, corresponding to the five symbols σ, α, β, η, and τ in Table 1. The flux lattices for them are illustrated in Fig. 1c-h, and are elucidated below. The technical details for how these flux lattices represent PSAs can be found in Methods.
(i) The first class refers to σ = L a L b L À1 a L À1 b for the translation subgroup. σ = ± 1 corresponds respectively to flux 0 and π through the plaquette spanned by L a and L b , as illustrated in Fig. 1c.
(ii) The second class concerns cohomology invariants α of symmorphic rotational symmetries. That is, α = R n for an n-fold rotation R, where α = ± 1 corresponds to flux 0 or π through the plaquette invariant under the rotation, as in Fig. 1d.
(iii) The third class corresponds to the square of a mirror reflection M, i.e., β = M 2 . It turns out that M 2 = − 1 cannot be realized on lattices with only nearest neighbor hoppings within one layer. We propose to realize it either by second neighbor hopping as in Ref. 28 or on a bilayer lattice with π flux through the interlayer plaquettes, as shown in Fig. 1e.
(iv) The fourth class (η invariants) is on relations between translations and reflections. For example, η = M x L y M À1 x L À1 y , and η = ± 1 corresponds to flux 0 or π through the plaquette in Fig. 1f x L À1 y = À 1, we may design the flux pattern as in Fig. 1g.
(v) The fifth class consists of invariants τ that extend the algebraic relations between translations and glide reflections, e.g., τ = g x L x g À1 x L x . As illustrated in Fig. 1h, τ = ± 1 respectively corresponds to flux 0 or π through the area spanned by the translation and glide reflection, which is half of the plaquette spanned by unit translations. It appears only for Pg group in Table 1.
With the above building blocks, we can systematically translate the cohomology invariants into fluxed lattices and obtain models realizing each of PSAs in Table 1. In this way, we have constructed a "canonical" lattice model for each wallpaper group G, in the sense that all PSAs for G can be realized in this single model, by simply varying the 0/π flux distribution in the lattice. In Methods, we categorize the 17 wallpaper groups into five classes to briefly introduce how the canonical lattice models are constructed.
As an example, consider P2 group with four α invariants. The algebraic relation for each α i in Table 1 corresponds to a twofold rotation center in the unit cell, as illustrated in Fig. 2a. Under lattice translation, each α i is associated with a class of translation-related rotation centers, which are distinguished by four colors in Fig. 2a. Then, the canonical model can be constructed with each plaquette hosting a unique rotation center (see Fig. 2a), corresponding to the dual lattice of the lattice of rotation centers. Each α i = ± 1 can then be realized by inserting flux 0 or π into the corresponding class of plaquettes.
The canonical models for all 17 wallpaper groups are illustrated in Fig. 3, and are explicitly constructed in Supplementary Note 3. For each wallpaper group in Fig. 3, the cohomology invariants correspond to independent fluxes in the lattice model, and we distinguish the fluxes by different colors. This is consistent with the number 2 n of PSAs, with n the number of colors in each lattice model.

Physical consequences of projective symmetry algebras
Our revealed PSAs can lead to a wealth of new physics, beyond conventional systems based on ordinary representations. Below, we present three remarkable consequences for demonstration.
(1) Shift of high-symmetry points. In ordinary band structures, high-symmetry points are located either at the center (Γ point) or on the boundary of Brillouin zone (BZ) 2 . In contrast, with PSAs, the highsymmetry points are redistributed, and they can be at non-central points in the interior of BZ.
For instance, continue with the example of P2 group. Let's consider the PSA with α 1 = α 2 = 1 and α 3 = α 4 = − 1 (see the canonical model realization in Fig. 2b). Clearly, in this case the two translations L a and L b commute as usual, and therefore the BZ is unchanged. However, since b by the fluxes, the R-invariant momenta are transformed from (0, 0), (0, π), (π, 0) and (π, π) to The first column 'G' lists the names of 17 wallpaper groups, and the second column 'H 2 ðG, Z 2 Þ' presents the corresponding second cohomology groups. The chosen generators of each wallpaper group are given in the third column. The PSAs and the cohomology invariants are presented in the fourth column. L, R, M and g denote translation, rotation, mirror and glide reflection, respectively. σ, α, β, η and τ are the five classes of cohomology invariants valued in Z 2 = f ± 1g. The last column 'N G ' lists the number of equivalence classes of PSAs for each wallpaper group.
(2) Enforced nontrivial Zak phase. While ordinary crystal symmetries may protect topological structures of energy bands, we discover that some PSAs can even enforce nontrivial topological structures. That is, once the PSA is realized, certain topological invariant is guaranteed to be nontrivial.
Here, we give one example of this fascinating phenomena, again using the P2 group. Let us consider the PSAs with R 2 = ðL a RÞ 2 = α and ðL b RÞ 2 = ðL a L b RÞ 2 = À α, which can be realized by the canonical model with the flux configuration in Fig. 2b. The PSAs lead to RL b R À1 = À L À1 b for both α = ± 1. From this relation, one can show that the anti-unitary operator RT will shift momentum k to k + G b /2 with G b the reciprocal translation vector corresponding to L b (see Methods). Now, consider the effect of RT on a single energy band with eigenstates ψ k . Recall that spacetime inversion symmetry can quantize the Berry phase, also known as the Zak phase, along any periodic path in the BZ to be either 0 or π 29 . In contrast, here, RT with (RT) 2 = α exerts a stronger constraint on the Zak phase, i.e., it completely determines the Zak phase θ b along any G b -periodic path as due to the nontrivial action of RT discussed above (see Methods). This result means: if α = − 1, the Zak phase is enforced to be nontrivial. This is confirmed by the concrete model in Fig. 2b. This model has four isolated bands, and each band is enforced to have a nontrivial Zak phase π along k b . Hence, there must be topological edge states within the first and the third energy gaps, as shown in Fig. 2d.
(3) Eightfold degenerate nodal point. Highly degenerate nodal points protected by crystal symmetries have been a hot topic. Without including the twofold degeneracy of spin-1/2, the highest degeneracy protected by wallpaper groups is fourfold 30 . Here, we find that PSA can achieve a degeneracy of eightfold, beyond any ordinary representations. This is exemplified by the PSA of P3m1 with L a L b L À1 a L À1 b = À 1 and M 2 = − 1 (see Table 1). The canonical model is illustrated in Fig. 2e and f. Since L a does not commute with L b , we choose L 2 a and L 2 b to generate an invariant subgroup of P3m1, and the BZ is specified by L 2 a = À e ikÁe a and L 2 b = À e ikÁe b under the Fourier transform, with e a,b being the translation vectors of L 2 a,b . At high-symmetry point Γ, the little co-group is given by Z 2 2 D 3 × Z T 2 , where Z 2 2 are generated by L a,b . This little cogroup is projectively represented with factors inherited from that of P3m1. We find that it has two 4D irreducible representations and one 8D irreducible representation. The latter gives the eightfold nodal point, which is indeed confirmed via a concrete model as illustrated in Fig. 2g and h.

Discussion
In conclusion, we have established a unified theory for T-invariant crystals. Particularly, we classified all PSAs of wallpaper groups, developed a general construction method, presented canonical models to realize each PSA, and revealed remarkable physical consequences. The theory can be directly extended to 3D space groups. Our work provides a solid foundation for the study of artificial crystals and opens the door to a wealth of new physics beyond the current paradigm based ordinary symmetry representations. Notably, although our focus here is on spinless systems (which most artificial crystals belong to), the generalization to spinful systems is straightforward. This is because in the presence of T-invariance, it is always sufficient to consider Z 2 -valued factor systems, as stressed above. Then, in addition to the phases arising from fluxes, one only needs to take care of reflections and rotations of spin-1/2 by 2π, which lead to the phase − 1. Hence, all the cohomology invariants in classes (ii) and (iii) are reversed. Formally, we may just replace each α and β by (−1) 2s α and (−1) 2s β, respectively, with s = 0 and 1/2 for spinless and spin-1/2 cases.
Finally, we note that our theory of PSAs is based on two fundamental principles of physics: (a) Physical systems are classified by symmetries (Landau's paradigm); and (b) Symmetries are projectively represented in a physical system (Wigner's principle). Hence, the PSAs derived here are general and classify all T-invariant crystal systems, including not only artificial crystals, but also real materials, strongly correlated spin systems, and beyond.

Projective representations with time-reversal symmetry
In the main text, we emphasized that with T symmetry, the phase factors of space group symmetries can be constrained to be valued in Z 2 . Here, we present a proof for this proposition.
Let us enlarge the space group G by including T with T 2 = 1. Then, each group element can be written as gT a with g ∈ G and a = 0, 1. Suppose that under the projective representation ρ, the phase factor λ arises through ρðg 1 T a 1 Þρðg 2 T a 2 Þ = λðg 1 T a 1 ,g 2 T a 2 Þρðg 1 g 2 T a 1 + a 2 Þ: ð3Þ We shall prove that by appropriately modifying the phase of each operator ρ(gT a ), we can always transform λðg 1 T a 1 ,g 2 T a 2 Þ into the form, where ν(g 1 , g 2 ), ωðT a 1 ,T a 2 Þ 2 Z 2 .

Projective symmetry algebras and gauge fluxes
Let us consider a set of lattice sites and hopping amplitudes among them, which give rise to a tight-binding model, Here, a y i and a j are the particle creation and annihilation operators at sites i and j, respectively. H ij represents the hopping amplitudes t ij from site j to i if i ≠ j and the onsite energy ϵ i at site i if i = j. H is a Hermitian matrix and called the one-particle Hamiltonian of the tightbinding model.
Each hopping amplitude t ij may have a phase e iϕ ij (such that t ij = |t ij |e iϕ ij ), which is called the gauge connection of the lattice model. Particularly, here we consider the Z 2 gauge connections with ϕ ij ∈ {0, π}. For each closed loop C formed by successive hoppings, one can compute the product W C of the phases of all the hopping amplitudes involved. W C is called the Wilson loop operator of the loop C, and the gauge flux Φ C through the loop C is given by W C = e ÀiΦ C . For the Z 2 gauge field, we have W C ∈ { ± 1} and Φ C = {0, π}.
For each site i, we may change the phase of a y i for each i by an arbitrary e iθ i . Particularly, θ i is valued in {0, π} for the Z 2 gauge field considered here. Accordingly, the hopping amplitudes are transformed as t ij 7 !e iθ i t ij e Àiθ j , which is called a gauge transformation. An immediate result is that W C = e ÀiΦ C is invariant under any gauge transformation. This can be seen from that the ending site of a hopping is the starting site of the next hopping in a loop C, and therefore all phase changes involved are cancelled out. To summarize, the gauge fluxes are gauge-invariant quantities, whereas the gauge connections are not.
Only gauge-invariant quantities are physical. In the current case, the gauge flux configuration completely determines the physics of the model. Hence, a spatial transformation R that leaves the crystal and the gauge flux configuration invariant is regarded as a symmetry of the system. However, R does not necessarily preserve the gaugeconnection configuration A. After the action of R, A is generally changed to another one A 0 . Since the two gauge-connection configurations A and A 0 describe the same flux configuration, they are related by a gauge transformation G R . On the lattice, R is represented by a matrix indexed by lattice sites, which we still denote by R. The gauge transformation G R is a diagonal matrix with ½G R ii = e iθ i R , i.e., with the phase assigned to the ith site. Then, the physical symmetry operator in this case should be the combination That is, after the spatial transformation R, the gauge transformation G R is needed to recover the original gauge connection configuration A.
Notably, it is R = G R R that commutes with the Hamiltonian H, i.e.,

½R, H = 0: ð15Þ
The commutation relation is equivalent to the requirement, where G R ðiÞ = e iθ i R , namely the phase assigned to site i, and R(i) is the site transformed from i by R.
Then, we consider the successive action of two spatial symmetries, R 1 = G R 1 R 1 and R 2 = G R 2 R 2 . There are two natural operators to implement the action, namely, G R 1 R 1 G R 2 R 2 and G R 12 R 12 with R 12 = R 1 R 2 .
is a diagonal matrix with ith diagonal entry being G R 1 ðiÞG R 2 ðR À1 1 ðiÞÞ=G R 1 R 2 ðiÞ, and therefore represents a gauge transformation. It is clear that Δ G (R 1 , R 2 ) commutes with all possible symmetry-preserving Hamiltonians. Particularly, let us presume the usual case that H is a connected lattice model, i.e., any two sites are connected by hoppings. The presumption sufficiently leads to the fact that Δ G (R 1 , R 2 ) is proportional to the identity matrix, namely ½Δ G ðR 1 ,R 2 Þ ij = νðR 1 ,R 2 Þδ ij with νðR 1 ,R 2 Þ 2 Z 2 & Uð1Þ, i.e., the physical symmetry operators satisfy the PSA, If ν and ν 0 are related by transforming R to R 0 = χðRÞR with χ(R) ∈ U(1) or Z 2 , the two PSAs belong to the same cohomology class. It must be noted that the cohomology class of such a PSA is independent of the choice of gauge connections, and is solely determined by the flux configuration.

Realization of cohomology invariants
Based on the general discussions in the last section, we now show the flux lattices in Fig. 1c-h can realize the five classes of cohomology invariants, respectively. (i) Let us start with the cohomology invariant σ = L a L b L À1 a L À1 b . Since the algebraic relation is equivalent to for any lattice site i. For the lattice model in Fig. 1c, we have from (16) the relations which implies Here, 1, 2, 3, 4 label the four sites in Fig. 1c, which are counted counterclockwise from the lower left corner. The flux through the rectangle satisfies e ÀiΦ = e iϕ 12 e iϕ 23 e iϕ 34 e iϕ 41 = e iϕ 12 e iϕ 14 e iϕ 21 e iϕ 41 G a ð2ÞG * Thus, in the presence of flux Φ, L a and L b satisfy L a L b L À1 a L À1 b = e iΦ . This argument can be generalized to other lattices. See Supplementary Figs. 28 and 37.
(ii) For a cohomology invariant α = R n 2π=n , n must be even. Here, we have added the subscript 2π/n for R to specify the rotation angle. When n is even, rotating n/2 times is a twofold rotation R n=2 2π=n = R π . In general, R n=2 2π=n = ξR π with ξ 2 Z 2 . No matter whether ξ = ± 1, the cohomology invariant can always be expressed as Substituting R π = G π R π into the identity above, we see the cohomology invariant can be realized by G π ðiÞG π ðR π ðiÞÞ = α ð24Þ for any site i. Let us label the vertices of the plaquette in Fig. 1d by i = 1, 2, 3, ⋯ 2l with n = 2l. Then, R π (i) = i + l. From (16), the hopping amplitudes satisfy t i,i + 1 = t i + l,i + l + 1 G π ðiÞG * π ði + 1Þ: Then, the flux is found to be From (24), we conclude that e iΦ = α. Note that all phases are restricted in Z 2 = f ± 1g.
(iii) For M 2 = β, the 2D mirror reflection is interpreted as the twofold rotation through an axis parallel to the 2D plane. Then, it is clear from (ii) that e iΦ = β with Φ the flux through each vertical plaquettes in Fig. 1e.
(iv) The cohomology invariant η = M x L y M À1 x L À1 y is translated as Hence, we have the identity, Label the four white sites in Fig. 1f by 1, 2, 3 and 4, which are counted counterclockwise from the lower left site. From (16), we have the identities, Then, the flux is computed as e ÀiΦ = e iϕ 12 e iϕ 23 e iϕ 34 e iϕ 41 = e iϕ 12 e iϕ 14 e iϕ 21 e iϕ 41 G m ð2ÞG * m ð3ÞG * y ð4ÞG y ð3Þ = G * m ð3ÞG * y ðMð3ÞÞG m ðL À1 y ð3ÞÞG y ð3Þ = η * : which is equivalent to For the lattice in Fig. 1h, the hopping amplitudes satisfy referring to (16). Here, 1, 2, 3, 4 label the four sites in Fig. 1h, which are counted counterclockwise from the lower left one. Then, the flux through the rectangular plaquette is e ÀiΦ = e iϕ 12 e iϕ 23 e iϕ 34 e iϕ 41 = e iϕ 12 e iϕ 14 e iϕ 12 e iϕ 41 G x ð2ÞG * Note again that all phases above are either + 1 or − 1.

Construction of canonical models
In the main text, we have elucidated how to realize cohomology invariants σ, α, β, η and τ in Table 1 by lattice flux configurations. Here, we briefly introduce the general procedure for constructing the canonical models for all 17 wallpaper groups, which realize all 458 PSAs. It must be noted that the purpose of these models is to demonstrate the physical realization of all PSAs, so they are made as simple as possible and contain only nearest neighbor hoppings. One can certainly write down more complex models with more complicated lattice structures and hopping processes for a given PSA, just like what one typically does when constructing models based on ordinary representations of space groups. The illustration for all the canonical models is given in Fig. 3, and the full details for the model construction can be found in the Supplemental Note 3. In this process, we find it useful to categorize the 17 wallpaper groups into five classes.
(a) Groups P1, P3, and Pg are quite simple, since each of them has only one cohomology invariant σ or τ (see Table 1). It is straightforward to design lattice models with the flux patterns as introduced in (i) or (iv).
(b) For groups P2, Pgg, P4, and P6, all cohomology invariants are of type α in (ii), i.e., each α i = R n for some n-fold rotation R through a rotation center in the unit cell. Under lattice translations, each rotation center gives a lattice of rotation centers. Accordingly, each α i is associated with such a lattice, and different α i^' s correspond to different lattices. This has been illustrated with our example P2 in the main text. Then, the canonical model is constructed as the dual lattice for the lattice of rotation centers. This means each plaquette in the canonical model hosts a unique rotation center; conversely, each rotation center is the center of a plaquette preserving the rotation symmetry. Then, each α i = ± 1 is realized by inserting flux 0 or π into the corresponding plaquettes.
(c) For groups Cm, P3m1, and P31m, each has two cohomology invariants σ and β. We first construct a one-layer lattice model realizing σ as described in (i). Then, we double the one-layer model into a twolayer model, and introduce the nearest-neighbor interlayer hopping amplitudes to realize β as given in (iii) or in Fig. 1e.
(d) Each of groups Pmm, Cmm, P4m, P4g, and P6m has the two types of cohomology invariants α and β in (ii) and (iii), respectively. Here, following (b), we first construct a one-layer model to accommodate all α-invariants. Then, we double the one-layer model into a two-layer model, and appropriately insert fluxes for interlayer plaquettes to realize all β-invariants. Note that according to Fig. 1e, the vertical mirror planes should cross lattice bonds rather than lattice sites, which can always be satisfied.
(e) The remaining two groups are Pm and Pmg, both having ηand β-invariants. Therefore, we refer to (iv) and Fig. 1g for constructing two-layer models for them. Since Pmg also has two α-invariants, we first construct the one-layer model according to the α-invariants following (b), and then double it into a two-layer model to accommodate the βand η-invariants following Fig. 1g.

Shift of high-symmetry points
We derive the shift of high-symmetry points in Fig. 2c. From the PSAs for P2, it is straightforward to derive that Alternatively, these relations can be derived from the configuration of the canonical model in Fig. 2a. There, the fluxes through the plaquettes colored in red, green, purple, blue correspond to cohomology invarants α 1 , α 2 , α 3 , α 4 , respectively. If α i = 1(α i = − 1), the corresponding flux is 0 (π). When α 1 = α 2 = α and α 3 = α 4 = − α, plaquettes in each row have the same flux, and the flux values alternate across the rows (see Fig. 2b). Accordingly, the third equation above gives RL b R À1 = À L À1 b . In momentum space, L b is diagonalized as e ik⋅b . Then, Re ikÁb R À1 = e ÀiðkÀG b =2ÞÁb , with G b the reciprocal lattice vector for b. Hence, we see that k is transformed to − k + G b /2 under R, i.e., Thus, the R-invariant momenta are shifted to ± G b /4 and G a /2 ± G b /4, as shown in Fig. 2c.

Enforced topology by projective symmetry algebras
Here, we provide the details for the nontrivial Zak phase enforced by RT symmetry discussed in the main text. The RT symmetry puts the following constraint for a Hamiltonian: where U RT is a unitary operator determined by RT. Suppose that HðkÞ has a single band jψ k i over the G b period from k = 0 to G b . The action of RT on jψ k i will give a band eigenstate at k + G b /2, generally differing from jψ k + G b =2 i by a k dependent phase, i.e., Accordingly, RT relates the Berry connection A b ðkÞ = ψ k |i∂ k b ψ k at k and k + G b /2 as Because of this, the Zak phase θ b = H dk b A b ðkÞ over any G b -periodic path can be expressed as Moreover, the PSA relation (RT) 2 = α requires that U RT U * RT = α, which in turn leads to Thus, we arrive at θ b = i ln α mod 2π as claimed in the main text. For α = − 1, this Zak phase is guaranteed to be nontrivial.

Eightfold degenerate nodal point
For the group P3m1, when σ = α = − 1, the little co-algebra at the Γ point has two irreducible 4D representations and one irreducible 8D representation. The 8D irreducible representation can be expressed as Here, K denotes the complex conjugation, and U R = expðin 1 Á σ2π=3Þ, D R = expðÀiσ 2 2π=3Þ, with n 2 = ð1, À 1, À 1Þ= ffiffiffi 3 p and n 2 = ð0,1, À 1Þ= ffiffiffi 2 p . Following (c), we can construct the canonical model for P3m1. First, we construct the one-layer lattice that realizes the σ-invariant for translations as illustrated in Fig. 2d. Then, we double the one-layer model into the bilayer model as illustrated in Fig. 2e. In order to realize all cohomology invariants, we add flux at all regular hexagons and rectangles. The expression for this lattice model is given in the Supplementary Note 4.

Data availability
The data generated and analyzed during this study are available from the corresponding author upon request.